Symmetry conditions for non-reciprocal light propagation in magnetic crystals 
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Recent studies demonstrated the violation of reciprocity in optical processes in low-symmetry magnetic crys- 
tals. In these crystals the speed of light can be different for counter-propagating beams. Correspondingly, they 
can show strong directional anisotropies such as direction dependent absorption also called directional dichro- 
ism[S. Bordacs et al., Nat. Phys. 8, 734 (2012); M. Saito et al., J. Phys. Soc. Jpn. 77, 013705 (2008)]. Based 
on symmetry considerations, we identify the magnetic point groups of materials which can host such directional 
anisotropies and also provide a list of possible candidate materials to observe these phenomena. In most of these 
cases, the symmetry of the crystal allows directional anisotropy not only for optical processes but also for the 
propagation of beams of particles and scalar waves. We also predict new types of directional optical anisotropies 
- besides the optical magnetoelectric effect and the magnetochiral dichroism investigated so far - and specify 
the magnetic point groups of crystals where they can emerge. 



When a measured quantity of matter remains invariant un- 
der the interchange of the source and the detector, the corre- 
sponding symmetry of the studied system is termed as reci- 
procity. Reciprocity in physics has a challenging definition. 
While Landau used the term reciprocity as a synonym of time 
reversal symmetry 1 , later it was generalized by also consid- 
ering spatial symmetries of the studied system 2 -. A rigorous 
definition of reciprocity for scattering processes has been re- 
cently given by Deak and FiilorA 

Specific to optical phenomena, among reciprocity-violating 
media we consider those which can distinguish between 
counter-propagating light beams irrespective of the polariza- 
tion state of the light, i.e. we say a material shows directional 
anisotropy when such non-reciprocity is present for each po- 
larization of the light beam. As a counter-example, we re- 
call the well-known Faraday effect when the non-reciprocity 
holds only for circularly polarized light 4,5 , while counter- 
propagating linearly polarized beams have the same transmis- 
sion and reflection properties. In this case the reciprocity is 
violated for certain but not all polarizations, thus, in our term 
the corresponding media show no directional anisotropy. 

Although the presence of such reciprocity-violating direc- 
tional anisotropy in the properties of bulk matter may seem 
strange and counter-intuitive, it has been recently verified in 
a large variety of low-symmetry materials. In an absorb- 
ing medium this directional anisotropy results in different 
strengths of absorption for beams traveling along opposite di- 
rections which is called directional dichroism. Similar direc- 
tional optical phenomena - schematically shown in FigJT] — 
have already been observed for emission 6 ' 7 , absorption 8 and 
diffraction 9 of light in materials belonging to two different 
symmetry classes where the directional optical phenomena 
are manifested in the corresponding two configurations. 

The optical magnetoelectric effect (OME) was observed in 
systems with finite ferrotoroidal moment 10 - which can arise 
in simple cases as the cross product of the static magnetiza- 
tion and electric polarization of the material, T = P X M - 
for light propagating along or opposite to the toroidal mo- 
ment vector. Strong OME was first detected by Hopfield and 
Thomas in the absorption spectrum of excitons in the polar 




FIG. 1: (Color online) Schematic representation of direction- 
dependent light propagation in (a) transmission, (b) reflection, (c) 
emission and (d) Bragg-scattering. In case of (a) transmission and (b) 
reflection, forward and backward propagating beams distinguished 
by the medium are simultaneously shown. ( c) Excitations by higher 
energy photons (blue arrow) decay via emission characterized by 
direction-dependent emission rates (red arrows), (d) Equal intensity 
for beams diffracted symmetrically to the direct beam is broken by 
directional anisotropy. 



CdS with wurtzite structure . Instead of reversing the light 
propagation, the sign of the toroidal moment was switched to 
be parallel and antiparallel to the velocity of light. A remark- 
able portion of the absorption was found to be odd function of 
both the magnetization and the electric polarization of the ma- 
terial as expected for OME. Later the same phenomenon was 
found for the polar ferrimagnet GaFe03 in the X-ray absorp- 
tion of the core electron excitations 1 1 and for absorption due to 
intra-atomic d-d transitions over the visible-ultraviolet spec- 
tral rang o 12 i 13 . Enhancement of OME in diffraction was ob- 
served using a grating made of the same compound 13 . Meta- 
materials composed of two dimensional ferromagnetic islands 
with polar shape can also host OME in diffraction 14 . 
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Materials without spontaneous electric polarization and/or 
magnetization can show OME as well if a toroidal moment 
T is induced by crossed static electric and magnetic fields, 
as was observed for the electronic f-f transitions of Er 3+ in 
Ei"i.5Yi.5Al50i2 15 - Besides the usual case discussed above, 
spontaneous toroidal moment can also be present in anti- 
ferroelectric and/or antiferromagnetic compounds, whenever 
t - Y,i r i x m i i s finite 16 , while T can be vanishing. Here, 
m, and r, respectively stand for the z'-th magnetic moment and 
its position. For example the spin texture present in skyrmion 
lattice*^— can produce such toroidal moment. 

Besides its manifestation in optical processes, similar di- 
rectional anisotropy has been reported in diffusive electronic 
transport. The drain-source resistance of symmetric field ef- 
fect transistors, when applying crossed magnetic field and 
electric field (via a finite gate-channel potential), showed dif- 
ference for current flowing along or opposite to the toroidal 
moment 20 . Directional anisotropy was predicted recently for 
the propagation of spin-waves in planar magnonic crystals 
with metallic over- and dielectric underlay ers 21 and in an ul- 
trathin Fe film on a W(110) surface 2 ^. The latter was ex- 
perimentally observed by spin-polarized electron energy loss 
spectroscopy 23 . In both cases ferrotoroidal moment was intro- 
duced parallel to the wave propagation by crossed magnetic 
field or spontaneous magnetization and electric polarization 
arising from the layered structure. 

Except for the case of CdS, the previous examples for 
OME, detected in absorption/reflection/diffraction or in elec- 
tric resistance, gave small relative directional differences in 
the range of 10~ 5 - 1CT 2 . Recent experiments on magneto- 
electric multiferroic materials show a relative OME in the or- 
der of unity as it was demonstrated for intra-atomic d-d tran- 
sitions of CUB2O4 in the near-infrared spectral regional and 
in the electromagnon spectrum of Ba2CoGe207 in the THz 
photon frequency range 25 , which may bring the possibility of 
their optical applications e.g. as magnetically controlled di- 
rectional light switches. 

The second type of optical directional anisotropy effect, 
the magneto-chiral dichroism (MChD) appears as a differ- 
ence in the absorption or emission of chiral systems for 
light beams propagating parallel and antiparallel to an ex- 
ternal magnetic field. The existence of MChD was first ex- 
perimentally proved for a luminescent transition of the chi- 
ral Eu((±)tfc)3 complex 6 * 7 - where the intensity of the emit- 
ted light slightly differed for propagation along or oppo- 
site to the external magnetic field. The effect was also 
shown in the infrared absorption of a-NiS04 ■ 6H2O2 6 . and 
in the enantioselective photochemical dissociation of chiral 
complexes 27 . Similar phenomenon, difference in phase re- 
tardation of counter-propagating light beams was found in 
chiral organic liquids 28,29 . The emergence of MChD was 
also demonstrated in Bragg scattering 9 . The electrical ana- 
log of MChD is the difference in conductivity of chiral sys- 
tems for currents flowing parallel and antiparallel to the di- 
rection of an external magnetic field. Such electric magne- 
tochiral anisotropy was observed in helically distorted bis- 
muth stripes 30 and also in chiral carbon nanotubes 31 . In all 
these cases the directional anisotropy was rather weak corre- 



sponding to a relative difference of 1CT 8 - 10~ 2 in the studied 
quantities. Similar to the case of OME, multiferroic materials 
can show also a significant MChD effect in both their infrared 
d-d excitations 32 and spin-wave absorption spectral. 

Motivated by the recently observed plethora of directional 
anisotropy effects and their potential in future applications us- 
ing multiferroic compounds, we study the general symmetry 
conditions allowing directional anisotropy in optical phenom- 
ena. Our focus is on crystalline materials where we system- 
atically identify those crystallographic magnetic point groups 
which can support directional anisotropy effects. Besides the 
symmetry conditions where OME and MChD can be realized, 
we found additional magnetic symmetry groups which can 
host new types of directional anisotropy effects. These new 
directional effects are specific to crystalline media in contrast 
to OME and MChD which can also emerge in the optical prop- 
erties of gases, liquids and amorphous solids. 

When neglecting the weak interaction light-matter interac- 
tion possesses time reversal 34 (r) and full (9(3) spatial symme- 
try, in other words transformation of the whole "light source- 
sample-detector system" by r ® 0(3) symmetry operations 
must not change the interaction process. Therefore, only the 
crystallographic and magnetic structure of matter can give 
constraints to the symmetry of the measured optical quanti- 
ties. Using Neumann's principle these constrains can provide 
relations between the tensor elements describing the linear op- 
tical response such as transmittance and reflectivity matrices. 
These relations can be experimentally probed by changing the 
polarization or propagation direction of light. We focus on 
such cases when the equality of these tensor components for 
forward and backward propagation is broken by the low sym- 
metry of the material. Thus, every symmetry elements, which 
could convert any of the optical processes into their reversed 
pairs, need to be eliminated from the magnetic point group of 
such crystal. 

According to our definition, directional anisotropy is al- 
lowed for a given propagation direction in a material, if 
no symmetry operation of the crystal can transform a light 
beam (of any polarization) traveling along this direction into 
a counter-propagating one without changing its polarization 
state. This means that the path of the electric polarization vec- 
tor determines the same curve in the laboratory frame upon 
forward and backward propagation. As a necessary condition, 
the material must lack spatial inversion and time reversal sym- 
metries because they convert counter-propagating beams with 
the same linear polarization into each other. To systematically 
determine the conditions of directional anisotropy in a crystal 
we have to consider each element of the corresponding mag- 
netic point group. Let us take a beam with arbitrary elliptical 
polarization propagating along the z-axis. Since the electric 
displacement and magnetic flux density vectors (D and B, re- 
spectively) lie in the xy-plane we only need to treat those or- 
thogonal transformations which leave this plane invariant. In 
crystals these are the n-fold proper (ri) and improper (77) ro- 
tations around the z-axis for n = 1,2,3,4 and 6 and around 
an axis laying in the xy-plane for n — 1 and 2. All of these 
operations can be combined with the reversal of time which 
will be denoted by a prime (')■ In this paper an improper rota- 
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tion means a proper rotation combined with spatial inversion, 
giving 1 for the spatial inversion and 2 for a mirror plane sym- 
metry. 

As long as we consider the propagation of unpolarized 
light, the necessary and sufficient condition for the presence 
of a directional anisotropy effect along the z-axis is that no 
symmetry of the system can connect the opposite propagation 
directions, i.e. transform light k f f z vector to -k it z. Half 
of the transformations listed in the previous paragraph fulfill 
this condition, since we have the freedom to combine a spatial 
operation with time reversal, which converts a propagation- 
reversing symmetry to a propagation-conserving one and vice 
versa. Proper rotations (n) and improper rotations combined 
with time reversal («') around the z-axis meet this criterion for 
any n e { 1 , 2, 3, 4, 6}. Besides that 2 and 2' rotations around an 
axis perpendicular to the z-direction also conserve the propa- 
gation direction. In magnetic point groups containing only 
these propagation-conserving operations - together with any 
other transformations which do not leave xy-plane invariant - 
the propagation along the +z and -z directions can be differ- 
ent not only for unpolarized light but for other unpolarized or 
scalar waves (e.g. acoustic waves) and for unpolarized beams 
of particles. 

Next, we also take the polarization state of light into ac- 
count, because there can be some propagation-reversing oper- 
ations which also alter the polarization state of light (for each 
polarization) and thus allow directional optical anisotropy. 
However, this will only slightly affect the results obtained 
above. Among the rotations around the z-axis only n and n' 
reverse the propagation of light. Both transformations rotate 
the axes of the polarization ellipse by — . While n' transposes 
clockwise and counterclockwise circulation of the polariza- 
tion, n preserves the sense of circulation in the laboratory 
frame. Since the direction of propagation is reversed, n' will 
preserve and n will change the the sense of circulation in the 
frame of the light beam. For n — 1 and 2, n' connects every 
elliptical polarization state with its counter-propagating pair, 
and n does the same for linearly polarized states. Thus, both 
prohibit directional optical anisotropy (in our term) as at least 
for specific polarizations the propagation direction is reversed 
and the polarization in the light frame is kept intact. For n — 5 
the «-times repetition of n' and n results in 1' and 1, respec- 
tively, and for n = 6 the | -times repetition of n' and n gives 2' 

and 2, respectively, which all forbid directional anisotropy. In 
contrast, for n — 4 no power of 4 can convert light beams with 
elliptical polarization to their oppositely propagating pair. For 
4' the same is true except for the right- and left handed circu- 
larly polarized states which are transformed into their counter- 
propagating pair. 

Among the operations around an axis lying in the xy-plane 
only 1', 1, 2 and 2 can reverse the propagation direction and 
preserve the polarization state of the light beam at least for 
some specific (e.g. linear) polarizations. For n — 21 even 
integers the Z-times repetition of an n-fold rotation (denoted 
as («) ) gives 2 and (n ) results either in 2 or 2 depending 
on the parity of I. For n = 41, («') 2 ' = («) 2 ' = 2 and for 
n = 21 + 1, (n') n = 1' and (n)" = 1. Therefore, the presence 



of 21, 27,4/', 4/, (21 + 1)' or 21+ 1 rotations around an axis 
within the xy-plane prohibits directional anisotropies along 
the z-axis. Rotations around arbitrary axes, which are nei- 
ther parallel nor perpendicular to z, cannot generally connect 
opposite propagation directions along the z-axis. 

To conclude the previous analysis, almost all symmetry op- 
erations which reverse the direction of the light propagation 
retain at least one polarization state and hereby they cannot 
be present in the magnetic point groups of crystals exhibiting 
directional anisotropy. The only exception is 4 which allows 
directional anisotropy for any light polarization except for un- 
polarized light where the difference cancels out. 

Table|T]lists the crystallographic magnetic point groups ex- 
cept those which contain the time-reversal or spatial inversion 
symmetry. We indicate which symmetry groups support direc- 
tional anisotropy along certain high-symmetry axes. In some 
magnetic point groups the presence of directional anisotropy 
can be straightforwardly demonstrated by showing a static 
symmetry breaking vectorial quantity of the material which 
transforms under the symmetry elements of the group in the 
same way as the propagation vector of light. 

A static toroidal moment (T or t) transforms under all pos- 
sible symmetry operations identically to the wave vector since 
the latter is the cross product of the electric and magnetic com- 
ponents of light. Therefore, systems with finite toroidal mo- 
ment vector allow directional anisotropy for beams propagat- 
ing parallel or antiparallel to this vector and the corresponding 
effect is termed as the OME. The magnetic point groups com- 
patible with T - denoted by T in Table|T]- are subgroups of the 
symmetry group of a perpendicular magnetic and electric field 
vector, 2'mm' , where the electric field is parallel to 2', mag- 
netic field is perpendicular to the mirror plane m and T points 
perpendicular to the plane m' . In the notation of the magnetic 
point groups we follow the international convention, where m 
stands for mirror plane symmetry. The toroidal moment t can 
arise without the presence of electric polarization and mag- 
netization, thus, t is compatible with higher symmetry point 
groups than T. Indeed, t can exist the same point groups as 
a general time-reversal odd polar vector. In the case of crys- 
tals these are the subgroups of n/m'mm and n'2'm where t is 
parallel to the principal axis and n — 4 or 6. Since a sponta- 
neous, non-dissipative electric current 62 has the same symme- 
try properties as f, they can be present in the same 31 crys- 
tallographic magnetic point groups. These groups are labeled 
with t in Table U 

Chiral systems only possess proper rotational symmetries, 
which may be combined with the reversal of time. All of 
these operations transform the axial vector of the magnetic 
field (or spontaneous magnetization) and the polar wave vec- 
tor in the same way. Therefore, in chiral systems with a time- 
reversal odd axial vector order parameter such as magnetiza- 
tion light propagation along or opposite to this vector can- 
not be connected by symmetry. The corresponding directional 
anisotropy is the so-called MChD and is indicated by MC in 
Table U 

Additionally to the well-known OME and MChD, we found 
if a polar magnetic point group contains only symmetry ele- 
ments which are not combined with time-reversal then these 
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TABLE I: Crystallographic magnetic point groups and possible configurations for directional anisotropy (DA). We omitted the point groups 
which contain the 1' time reversal or 1 spatial inversion operation. After the international symbol of each point group the symmetry allowed 
magnetic structure (F - ferromagnetic, AF - antiferromagnetic) is indicated. The potential directional anisotropy configurations are listed for 
wave propagation parallel and perpendicular to the principal axis as || to PA and _L to PA, respectively. Here i) T and t, ii) MC, and iii) P mean 
the directional anisotropy effects emerging in i) ferrotoroidic materials for propagation parallel to the toroidal moment T and t, ii) in chiral 
systems for light traveling along the magnetic field, and iii) in point groups, where none of the group elements is combined with the time 
reversal for propagation along the ferroelectric polarization, respectively. X denotes cases where directional anisotropy can also be present but 
it cannot be characterized by a static vector quantity. Several examples of multiferroic crystals belonging to the given magnetic point group 
are given in the 5th and 10th columns, where B a stands for an external magnetic field pointing to the a crystallographic direction. 
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32 (AF) 
32' (F) 
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t 
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t 
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X 
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P t 
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A' lm' (AF) 
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- 
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42' m' (F) 
4mm (AF) 
A' mm' (AF) 
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X 






4' 2(4' 2' 2) (AF) 










m'3 (AF) 


X 






42' (42' 2') (F) 


MC, t 




NdsSi^ 4 . 
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43m (AF) 
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A' 1 m'm'm (AF) 










4 3m' (AF) 










A/m'mm (AF) 


t 








43(432) (AF) 










A/m'm'm' (AF) 










4'3(4'32) (AF) 
m'3m (AF) 


X 
X 
















m'3m' (AF) 









"FeGai£, MnSii2 B [m] ; RbFe(Mo0 4 )2^ 
i, Co 3 B 7 0i 3 Cl 3 i 
c BiTeI B [100 ^ 
<, FeGe 18 ,MnSi l<) B [na] 

e Cu 2 OSe0 3 B[iioA Co 3 Te0 9 4 ^;NdFe 3 (B0 3 ) 4 B [ow ^. 
'MnPS 3 46 ; Co 3 Te0 6 17 K < T < 21 K£ 
«Nb 2 Co 4 09^ 

''CuB 2 4 Bino] 24 ; CdS 8 , A1N, GaN, InN B [l0Q fiL; CaBaCc^OA 
GaFeOjil; Co 3 B 7 Oi 3 Bi^; KMnFeF^ 

'CuB 2 0<A Ca 2 CoSi 2 7 ^, C^OSeOjil B [10 oj; 

[Ru(bpy) 2 (ppy)][MnCr(ox)]£i 
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operations transform the vector of electric polarization or tem- 
perature gradient, both even in time reversal, in the same way 
as the wave vector being odd in time reversal. In this case the 
presence of a static polarization/temperature gradient makes 
the wave propagation along or opposite to it inequivalent. This 
potential directional anisotropy effect, denoted by P in Table 
H] has not been studied either theoretically or experimentally 
so far. 

Please note that whenever symmetry permits directional 
anisotropy of either MC- or P-type, t-type directional effect 
can also emerge. Nevertheless, the distinction between these 
cases is meaningful as the physical origin of these directional 
effects can be different. Besides these cases we found addi- 
tional situations when the symmetry condition for directional 
anisotropy cannot be related to any static vectorial quantity 
which transforms in the same way as the wave vector. These 
point groups are marked with an X in Table U 

Except for the magnetic point group 4, all the cases of di- 
rectional anisotropy listed in Table J] also have validity for 
other processes than light propagation. Symmetry elements of 
these groups cannot reverse wave propagation and, thus, allow 
directional anisotropy for propagation of electrons, neutrons, 
acoustic- and spin-waves, and in electric or heat conduction. 
In some other groups directional anisotropy is conceivable but 
only for polarized light, spin polarized neutron and electron 
beams or any transversal wave. These cases can be determined 
by systematical checking whether the propagation-reversing 
operations can preserve the polarization state. Specific to light 
propagation, we found that 4 is the only additional point group 
supporting directional optical anisotropy. 

Another approach to study directional optical anisotropy 
is the symmetry analysis of the optical response functions. 
Due to its simultaneous time-reversal and spatial inversion 
breaking nature, optical directional anisotropy can only be 
produced by the bilinear magneto-electric tensor when dipole 
approximation appropriately describes the light-matter inter- 
action. The non-vanishing elements of this tensor can be 
identified by the symmetry analysis of the given material. 



For the static magneto-electric tensor the result of such anal- 
ysis based on magnetic point groups can be found in the 
literature 63 . While the static magneto-electric tensor can only 
contain time-reversal odd elements, the dynamic one may 
have time reversal even part as well. Therefore, the sym- 
metry dictated form of the dynamic magneto-electric tensor 
needs to be included in the constitutive relations in order to 
solve the Maxwell's equations and check the presence of di- 
rectional anisotropy. This approach is specific to optics and 
its validity is restricted to the dipole approximation of light- 
matter interaction. Considering compounds with cubic mag- 
netic point groups, in dipole order no directional anisotropies 
can be present since the magneto-electric tensor is propor- 
tional to the unit tensor. However, taking into account higher 
order terms described e.g. by electric quadrupole tensor, the 
symmetry allowed directional anisotropy along the threefold 
axis shown in Table U can be reproduced. Nevertheless, we 
think that our approach provides a more general framework to 
address the question of directional anisotropies in crystalline 
solids. 

In conclusion, possible direction-dependent optical config- 
urations were determined using symmetry considerations. Be- 
sides the conventional OME and MChD a third potential di- 
rectional anisotropy effect was predicted when the wave prop- 
agation is parallel to the electric polarization in materials 
whose symmetry groups do not contain any symmetry opera- 
tion involving time reversal. With one exception, all the listed 
configurations allow non-reciprocity for any transport phe- 
nomenon, particle or wave propagation. We identified each 
crystallographic magnetic point group which can host direc- 
tional anisotropies providing a key for the systematic investi- 
gation of candidate materials. We also propose a large collec- 
tion of materials where these effects are expected to emerge. 

We thank G. Kriza, Y. Tokura, L. Deak, T. Fiilop and B. 
Dora for valuable discussions. This project was supported 
by Hungarian Research Funds OTKA PD75615, TAMOP- 
4.2.2.B-10/1-2010-0009. 
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